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Abstract 

In this paper, it is shown that Bermudan option pricing based on either the reduite (in a one¬ 
dimensional setting: piecewise harmonic interpolation) or cubature - is sensible from an economic 
vantage point: Any sequence of thus-computed prices for Bermudan options with increasing sets 
of exercise times is increasing. Furthermore, under certain regularity assumptions on the payoff 
function and provided the exercise times are equidistant of exercise mesh size h, it has a supremum 
which coincides with the least fixed point of the approximate pricing algorithm - this algorithm 
being perceived as a map that assigns to any real-valued function / (on the basket of underlyings) 
the approximate value of the European option of maturity h and payoff function /. 



1 Introduction 


The Bermudan pricing approach outlined in the first part of this paper is to iteratively construct 
a piecewise harmonic approximation to the function that assigns the value of a Bermudan option 
(with payoff function g, exercise mesh size t > 0 and maturity nt) to the vector of logarithmic 
start prices of the underlying assets. 

In the first step of this iteration, one will compute a piecewise harmonic approximation to the 
function that assigns the corresponding European option price of maturity t to the logarithmic 
asset prices at the penultimate time T—t where exercise is possible. The iteration step will consist 
in computing the expectation of this function after time t (under the assumption that the vector 
price process is a Markov process), discount, take the maximum with the payoff function g , and 
then perform a reduite-based interpolation (in the one-dimensional setting: an interpolation by 
interval-wise harmonic functions). (This procedure computes a piecewise harmonic approximation 
to the so-called Snell envelope (2j from the theory of optimal stopping Eiffl) 

The notion of harmonicity will be derived from the log-price process X, assumed to be a 
time-continuous Markov process with a path-continuous modification, i e a diffusion - which 
then makes the infinitesimal generator of X a second-order elliptic differential operator and gives 
rise to a notion of harmonicity. 

The principal advantage of such an algorithm based on piecewise harmonic interpolation is 
that it is intrinsically finite-dimensional, since the space of harmonic functions corresponding to 
X will be two-dimensional. However, one needs to make sure that this algorithm will make sense 
economically. In particular one will have to require that the iteration step be monotone, since 
adding one more possible exercise time of course increases the value of the option price. This 
is not difficult to see and will be one of the first results of this paper. Also, once monotonicity 
of the iteration step has been established, the question of whether there exists a supremum to 
this sequence of approximate non-perpetual option prices arises, since this could then be treated 
as an approximation to the price of the corresponding perpetual option of exercise mesh size t. 
Under certain assumptions on the payoff function, this will be proven as well; furthermore we will 
be able to characterise it as the minimal fixed point of the iteration step previously referred to. 

First, we will discuss these questions in the one-dimensional setting - very little knowledge 
of potential theory has to be assumed for the proofs in that section. Second, we shall generalise 
that approach to higher dimensions; this will entail a few technical subtleties. At the end of both 
Sections we will apply these results to financially relevant settings. 

Finally, in a last Section we shall discuss related questions for Bermudan pricing based on 
cubature and establish a linear convergence rate that corresponds to the discount factor. 

2 Some terminology for Bermudan option pricing algo¬ 
rithms 

Let C° (R d , [0, +oo)) and L° (R d , [0, +oo)) , as usual, denote the spaces of nonnegative contin¬ 
uous functions defined on R d , and of equivalence classes of nonnegative measurable functions 
defined on R d , respectively. 

The observed price for any financial derivative increases when adding one more possible exercise 
time; also, the perpetual limit often has to be finite (in particular for puts, since then any option 
price is bounded by the strike price). 

Approximative Bermudan pricing algorithms also should reflect these features of actual ob¬ 
served prices, therefore we introduce the following manner of speaking to distinguish economically 
sensible from less sensible algorithms. 

Definition 1. A map D : L° (R d , [0, +oo)) —> L° (R d , [0, +oo)) is said to be a sound iterative 
Bermudan option pricing algorithm (for short, a sound algorithm,) for a payoff function 
g £ C° (R d , [0, +oo)) if and only if Df > g for all f € L° (R d , [0, +oo)) and the map D is 
pointwise monotone, that is 

V/o,/t G L° (K d ,[0,+oo)) 

(Vtel d fo{x)<h{x)=>\/xeR d Dfo(x) < Dh(x)) . 
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A sound iterative Bermudan option pricing algorithm D is said to have a perpetual limit if 
and only if 

sup D° n g G L° (R d , [0, +oo)) 

n£No 

(rather than this supremum being allowed to equal +oo on a subset of positive Lebesgue 
measure of its range). In that very case, the function in the last line is simply referred to 
as the perpetual limit of the algorithm. Finally, D is said to converge linearly in L°° to the 
perpetual limit if and only if there exists a c G (0,1) such that 

V„ 6 N || (D-+‘ - D") s|| L . ( „„ < c ||(D“ - B- 1 ) <?|L_ (RJ , R) , 

D n being shorthand for D° n for all n G No- 

Remark 1. The elements of L° ()R d ,[ 0,+oo)) should be conceived of as assigning the value 
- that is, the expected payoff - of an option to the vector of logarithmic start prices of the 
components of the basket (at least on the complement of a Lebesgue null set). 

Remark 2. The monotonicity condition imposed on sound iterative Bermudan pricing al¬ 
gorithms entail that the sequence of functions (D n g) nGNo = (D on g ) neNq is always pointwise 
increasing. Thus, this sequence has a limit: 

g < sup D n g = lim D n g. 

nGN n—>oo 

The infimum of all H-fixed points is always an upper bound for the perpetual limit: 

Lemma 1. Let D be a sound iterative Bermudan pricing algorithm for g with a perpetual 
limit u. Then the function u is smaller than any fixed point of D; moreover, Du > u. 

Proof. Any fixed point h of D is in the image of D and therefore, due to our assumptions 
on sound algorithms, pointwise greater than or equal to g. Now, as D (and thus D n ) is 
pointwise monotone, 

Vn G N 0 h = D n h > D n g , 

therefore 

h = sup D m h > sup D m g, 

m m 

where the right hand side is just the perpetual limit. Hence, any fixed point of D is greater 
than or equal to the perpetual limit. Furthermore, observe that due to the pointwise mono¬ 
tonicity of D , 

D (sup D n g S j > D m+1 g , 

therefore for all m G No, 

D (sup D n g) > sup D m+1 g = sup D m g. 

\ n J m m 

□ 

Later on, it will turn out that if D is based on piecewise harmonic interpolation or the reduite, 
the perpetual limit is, in fact the minimal fixed point (cf Theorem |T| and Lemma 0for piecewise 
harmonic interpolation, Theorem |2] for reduite-based approximation, and Theorem 0for a result 
on a map D which is based on cubature). 

3 Piecewise harmonic Bermudan option pricing for op¬ 
tions on one asset 

3.1 Definitions and facts from classical potential theory 

Let (P*)t>o be a Markovian semigroup associated to a one-dimensional diffusion and let L be the 
infinitesimal generator of L. Then L will be a second-order differential operator (cf eg Revuz and 
Yor 0). 
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If one assumed, for the sake of an example, that L is a second-order differential operator even 
with constant coefficients, that is 


L:f~pf' + \f", 

then P will merely be the semigroup of finite-dimensional distributions of a multiple of Brownian 
motion with linear drift. 

Let now U C I be open and nonempty. A twice differentiable function / : U —> R is 
called harmonic on U if and only if Lf = 0 on U. A continuous function g : U —> R is 
said to be subharmonic on U if and only if Lg is defined as an element of R U {± 00 } and 
Lg{x) = lirn t |o Pt9 G)-9i x ) > 0 for all x G U. A function / : U —*■ R will be called superharmonic 
on U if and only if —/ is subharmonic on U. 

A function that is harmonic (subharmonic, superharmonic) on R will simply be called harmonic 
(. subharmonic, superharmonic). 

Also, if h is both superharmonic and subharmonic, then one will have h > Pth > h for all 
t > 0, making h harmonic. 

The space of all harmonic functions on U is a vector space. As a classical result (cf eg Protter 
and Weinberger 0), one has that all subharmonic functions obey the following 

Maximum Principle: If g is subharmonic on U C R and U is compact, then the maximum of 
g will be attained on dU . 

Also, the set of subharmonic functions has the following closure properties: 

Subharmonic functions form a cone, closed under V and P: If f,g : R —>• R are subharmonic 
functions, then so are fVg, Ptf for all t> 0, and af + fig for all a, (3 > 0. 

See textbooks like Revuz and Yor pj, Meyer (33, Port and Stone [7} and the classical references 
therein. 

3.2 Harmonic interpolation 

Lemma 2 . Given two support abscissas ao ^ a\ and ordinates Co, c±, there is a unique 
harmonic interpolation, that is a harmonic function h : R —► R such that hfap = Ci for 
i G {0,1}. Thus, there cannot be more than two linearly independent harmonic functions, 
and two harmonic functions coincide on all of R, once they agree in two points. 

Assume now, for the sake of simplicity, that L is a second-order differential operator with 
constant coefficients, that is 

L:f~0f + \f". 

Then, if P ^ 0 , all harmonic functions have the shape 70 + 7 i • exp (-2/3-) for some real 
constants 70 , 71 , and if P = 0 , the vector space of all harmonic functions is just the space of 
all affine functions. 

Proof. The existence and uniqueness assertion for harmonic interpolation is merely the state¬ 
ment of the unique solvability of the Dirichlet problem for L. 

If there were more than two linearly independent harmonic functions on (ao,ai), the 
operation of harmonic interpolation would not be unique. 

Finally, constant functions are of course harmonic. Also, in case P = 0 any linear function 
must be harmonic as well, whence id« and 1 are already to linearly independent functions 
for the case P = 0. For the case of non-zero P, we remark that 

L exp (—2 P-) = —2/3 exp (-2/3-) • /3 + 4/3 2 exp (—2/?•) • i = 0 

and conclude that 1 and exp (—2/3-) are two (obviously linearly independent) harmonic func¬ 
tions. 

□ 

Lemma 3. Let I = (a, b ) be an open nonempty interval, and let f : I —► R be subharmonic 
on U such that f(a) = f(b) = 0. If there exists an xq E I with f(x 0 ) = 0, then f = 0 on I. 
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Proof. Due to the Maximum Principle, / < 0. Now suppose f(x o) = 0 for some xq G (a, b). 
If neither / = 0 on (a, xq) nor / = 0 on (xo, b), then we would have a yo G (a, xq) and a yi G 
(x 0l b) such that f{y 0 ) < 0 and f(yi) < 0. Then max 9 ( y0iyi ) / < max [y0iyi ] / = }(x 0 ) = 0, a 
contradiction to the Maximum Principle. Hence, either / = 0 on (a, xq), or / = 0 on (xo,b), 
and without loss of generality, we shall from now on assume / = 0 on (a, xo). 

If there was now an Xi G (xo,b) with f(x i) < 0, pick any a\ G (a,x o) and let (by virtue 
of Lemma 0) hi ■ R -► R be the unique harmonic function such that hi(ai) = f(ai) = 0, 
h\{x\) = f(x i) < 0. Then / — hi is subharmonic on (ai,cci) and vanishes on the boundary 
of (ai,£i). Hence, by the Maximum Principle applied to / — hi, one has f — hi < 0 
on (oi,a;i). However, / = 0 on (ai,iro) C (a,Xo) as established previously, and, by the 
Maximum Principle applied to hi, also hi < 0 on (ai,xi) D (ai,a;o). Thus 

0 < —hi = f — hi < 0 on (ai, xq), 

which means hi = 0 on (ai, Xq). But since (oi, Xo) is a nonempty open interval (and therefore 
does not have less than two elements) and hi is harmonic, this can only be true if hi is 
identically zero on all of R, in contradiction to hi(xi) < 0. 

□ 

If there is only a first order non-zero term, the space of harmonic functions will just coincide 
with the space of constant functions. 

Lemma 4. Let ao < ■ ■ ■ < a m G R. Consider a function f : R —> R that is subharmonic 
on ( ao,a m ) and let, by virtue of Lemma\^ 1(f) be the unique interpolation to f such that 
1(f) = f on {do,..., a m } and for all i < m, 1(f) coincides with some ha7'monic function on 
(a u a i+ 1 ), denoted by fa. 

Then the following assertions hold: 

1■ 1(f) > f on [a 0 ,a m \. 

2. 1(f) = max{/o,..., / m _i} on (ao,a m ), thus making 1(f) subha7'monic on (ao,a m ). 

3- 1(f) < f on (— oo, a 0 ) U (a m ,+oo). 

Corollary 1. Piecewise harmonic interpolation with respect to a set of support abscissas 
{do, ... ,a m } preserves subha7'monicity on (ao,a m ) in the following sense: If, in the notation 
of Lemma ^| / is subha7'monic on (ao,a m ), then so is 1(f). 

Proof of Lemma EJ 1. We apply for every i < m the Maximum Principle to the subhar¬ 
monic function / — /,; on (aj,ai+i) (which vanishes on the boundary of (aj,dj +1 )) to 
get 0 > f - fi — f - 1(f) on [at, a i+ 1 ] for all i < in. This proves 1(f) > f on [a i? a i+ i]. 

2. Again from the Maximum Principle, we get that two different harmonic functions can 
only have one point in common. Hence, not only {fi = fi+i} 2 {oj+i} for all* < m— 1 
(a trivial consequence of the interpolation) but if /., ^ fi+i, then even 

{fi = fi+l} = { a i+l} ■ 

Consider any i < m. There are two possibilities: Either fi < f i+ i on (—oo,a,; + i) and 
fi > fi+i on (a i+ i, +oo) or the other way round fi > f i+ i on (-oo, a i+ 1 ) and fi < f i+ i 
on (cq+i, +oo). 

Suppose, for a contradiction, the former situation holds for some i < m: fi < fi + i on 
(— oo, di+i) and /* > f i+ i on (a i+ i, +oo). Then 1(f) would equal fiAf i+ i on [ ai,a i+2 ], 
which is superharmonic. Then, / — (/,; A fi+i) would be subharmonic on [ai, 2 ] and 
it would have three zeroes, in ai, dj+i and a,j+ 2 . By Lemma El this can only be true 
if / - (fi A f i+ 1 ) = 0 on all of [ai,a i+2 \. Thus, f = fi A f i+ 1 on [ ai,a i+2 ]■ Since / is 
subharmonic on [ai,ai + 2 ], so must then be fi A fi + 1 (which is already superharmonic 
on [aj,cq+ 2 ]) then, and therefore, fi A fi+i is harmonic on [ai,ai+ 2 ]. Now, fi and /i+i 
agree only on one point, viz. aij+i, inside (ai,ai + 2 ). Therefore, /,; A fi + \ has with both 
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fi and fi- )_i at least two points in common. Since any harmonic function is uniquely 
determined by a mere two points, we thus get 

fi = fi A fi -\-1 = fi+1 
in contradiction to our assumption. 

Therefore, fi > fi+i on (—oo,a.j+i) and fi < /i+i on (a,;+i,+oo) for all i < m. 

Thus, due to ao < • • • < a m , we gain 

fo > fi 
fi > h 


fm —2 ^ fm—1 

as well as 

fo < fi 

h<h 


fm —2 ^ fm—1 

Hence we arrive at 

fo > fi > ■ ■ ■ > fm—i on (a 0 , ai) 

as well as 

Vi G {1..., m - 1} (fo < ■ ■ ■ < fi: /*>•••> fm- 1 ) on (ai, a i+ 1 ) 
Therefore, 

Vi < m Vj fi> fj on (a i; a i+ i), 

so 

Vi < m fi= max {/ 0 ,..., f m -i} on (a*, a i+ 1 ) 
which proves 1(f) = max{/ 0 ,..., f m - 1 }- 

3. This part of the proof is similar to the proof of Lemma 0 The function / - 1(f) 
is subharmonic on (—oo,ai) and it has two zeroes in ao and a\. First consider the 
case where /, whilst being subharmonic, is not harmonic on (oo,ai), which ensures 
that the subharmonic function / — 1(f) - which vanishes on the boundary of (ao,ai) 
- is not harmonic. Then there is an aq £ (ao,ai) such that (f — 1(f)) (x i) < 0. If 
/ — 1(f) was now also negative for some Xq £ (—oo,ao), then the Maximum Principle 
would yield 0 > max{ l0ill ) (/ — 1(f )) = max[ X0!Xl ] (/ —1(f)), but on the other hand 
(/ —1(f)) (ao) = 0 and ao £ (xo,Xi), a contradiction. Hence in the case where / is 
not harmonic on (ao, ai), we already have / — 1(f) > 0 on (—oo, ao). 

Next suppose / is indeed harmonic on (ao,ai). Then, so will be / —1(f) and because 
/ — 1(f) = 0 on {ao,ai}, we will get / — 1(f) = 0 on all of [ao,ai]. Now pick 
some Xi £ (ao,ai) and assume, for a contradiction that (/ —1(f)) (xq) < 0 for some 
Xq £ (—oo, ao). According to LemmaEJ there exists a unique harmonic function h such 
that h( x 0 ) = (f -1(f)) (xo) < 0 and h(x i) = (f -1(f)) (aq) = 0. Then f-l(f)-h is 
subharmonic on (xq, X\) and vanishes on the boundary of (a’o, aq), hence f—l(f) — h < 0 
on [xoj^i]- But also, h is harmonic and nonpositive on the boundary of (xo,X\), thus 
—h will be nonnegative on all of (cco,aq). Therefore 

0 < —h = f — 1(f) — h < 0 on [ao, aq], 

therefore h = 0 on (ao,aq) ^ 0, yielding h = 0 everywhere. This contradicts h(x o) < 0 
and hence completes the proof for / — 1(f) > 0 on (—oo, ao). 

Analogously, one can prove the domination of 1(f) by / on (a m , +oo). 

□ 


on (—oo, ai) D (—oo,ao), 
on (—00, 02) D D (—00, ao), 

on (-00, a m _ 1) ■ ■ ■ D (-00, a 0 ) 

on (ai, +00) D ■■ ■ D (a m , 00), 
on (a 2 , +00) D ••• D (a m , 00), 

on (a m _i,+oo). 
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Lemma 5. Piecewise harmonic interpolation to a set of support absicssas {ao,...,a m } is 
monotone on [ao,a m ] in the sense that if f < g on [ao,a m ], then the piecewise harmonic 
interpolation X(f) of f will be dominated by the piecewise harmonic interpolation X(g) of g 
On [ao, O'm] • 

Proof. Consider any i < m. By deflation, X(f) = / < g = X(g) on {aj,ai+i}, hence 
1(f) —1(g) < 0 on {a;,ai+i}. But T(f) —1(g) is harmonic on (aj,ai+i), therefore the 
Maximum Principle yields 1(f) —1(g) < 0 on [a,, a 1+ i]. Since this holds for all i < m, we 
get 1(f) — 1(g) < 0 on all of [ao, a m \ □ 

Lemma 6. Let, as before, X : R[“o> a ™] —> R® denote the operator of piecewise harmonic in¬ 
terpolation with respect to the set of support abscissas {ao,..., a m }. Consider a subharmonic 
function f : R —> R and a harmonic function h : R —> R such that f < h on R. Then 
1(f) < h on R. 

Proof. From the previous Lemma 0 we already know that X(f)(x) < X(h)(x) holds for all 
x £ [ao, a m ]. However, X(h) = h , hence X(f) < h on [ao, a m ] and from LemmaQJ we conclude 
that h> f > X(f) on the intervals (—oo, ao) and (a m , +oo). □ 

3.3 A fixed point theorem 

Theorem 1. Let X : R[ a °>°™l R r again denote the operator of piecewise harmonic inter¬ 
polation with respect to the set of support abscissas (ao,..., a m }, and fix t > 0, r > 0. Let g 
and c be subharmonic functions, and let h be harmonic and nonnegative. Let c be, moreover, 
ha7'monic on each of the intervals [aj,ai+i] for i < m, and assume Ptc(x) is finite for all 
i£l. Suppose c < g on [ao, a m ] and c,g<h on R. Now define 

1C: f ^X (e~ rt P t (1(f) Vc)V?) [ [a 0 , a m \ 


as well as 


{ / : R —> R subharmonic, f > c on [ao,a m ], 

/ [ [ao, a m ] : V* £ {1,..., m — 2} / harmonic on (a*, Oj+i), 

/ harmonic on (—oo, ai), (a m _i, +oo), / <h 

Then 1C maps the convex and bounded subset Q ofC°[ao, a m ] continuously to itself. Moreover, 
due to Lemma [U Q is a subset of a finite-dimensional subspace of C°[ao,a m ] (this subspace 
being the space of all functions from [ao,a m ] that are harmonic on each of the intervals 
[ai,cii+ 1 ] for i < m. By Brouwer’s Fixed Point Theorem, 1C has got a fixed point in Q. 
Finally, 1C is a composition of monotone functions on [oo,a m ] and therefore monotone as 
well. 

Proof. We can divide the proof for 1C(Q) C Q into three parts: 

1. The cone of subharmonic functions is closed under V, under Pt, under multiplication 
by constants and under piecewise harmonic interpolation X (cf LemmaQJ, therefore the 
image of Q under K, can only consist of subharmonic functions. 

2. By the monotonicity of P t and X (Lemma 0 , combined with the equations T\h = h 
and X(h) = h, as well as c,g< h, we have for all / : R —> R such that / \ [ao, a m ] £ Q, 




X g = e~ 


Pth \Jg < h. 
=h> o 


But the left hand side of this last estimate is subharmonic (because the cone of sub¬ 
harmonic functions is closed under V, Pt and, as we have remarked in Corollary^ also 
X). Hence, due to LemmaEl we obtain ICf < h. 
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3. The lower bound follows again from the monotonicity of 1 (Lemma 0 that is), but this 
time only by exploiting c < g on [ao, a m ] and the specific shape of c: For all / : M —> R 
such that / I" [ao,a m ] E Q, one has 

ICf = !(■■■ Wg)> 1(g) > 1(c) = c. 


We get that Q is bounded by supu 0)O 1 h V |c| > 0 as a subset of C°[ao, a m ], and because Q 
is finite-dimensional, we may apply Schauder’s Theorem, as soon as we have established the 
continuity of K, with respect to the norm j| ■ ||q : / i—> max{|/(ao)|,..., \f(a m )\}. 

Note, for this sake, that from the Maximum Principle, we have ||/||q = max[ aoam ] |/| 
whenever / is harmonic on each of the intervals (dj,aj+i for i < m. Therefore, if (/ n )neN G 
Q N , / G Q and ||/n - f\\ Q —* 0 for n -> oo, then f n —> / uniformly on [a 0 ,a m ] as 
n —> oo, hence - since all f n and / are harmonic on (— oo, ai) and (a m _i, +oo) - by Harnack’s 
inequality (cf eg |HJ p 11 If]), fn —► / pointwise on all of R as n —> oo. But 

c(x) < f n (x) < h(x) = P t h(x) 

for all n G N and iel, therefore Lebesgue’s Dominated Convergence Theorem (applied to 
the measure A i—> Pt\A(%) for each x G K) may be employed to get Ptf n (x) —> Ptf(x) as 
n —> oo for all ifl Via 


< 

< 

< 


e~ rt P t (1(f) V c) V g - e~ rt P t (l(f n ) Vc)Vj 
e~ rt P t (1(f) V c) - e~ rt P t (l(f n ) V c)| 
e~ rt Pt (1(f) V c — l(f n ) V c) | 
e~ rt P t (1(f)-l(f n ))\ 

e~ rt Pt (f — fn) I ) 


we can now deduce that \ICf n — ICf ||q —> 0 as n —» oo. 


□ 


The existence of a minimal fixed point for JC can be proven constructively as well: 

Corollary 2 . Let us adopt the notation of the previous Theorem. Then the sequence 
(£”(<? v 0)) neNo is monotone on [ao,a m ], bounded and dominated by h. Therefore we have 
the existence of a limit on [ao,a m ] given by 

\/x G [ao, a m ] q(x) := lim K. n (g V 0)(a;) = sup lC n (g V 0)(a;). 

n—> °° neN 0 

This limit is an element of Q and therefore can be canonically extended to the whole o/K. By 
the continuity of 1C, q is a fixed point of 1C. On [ao,a m ], the convergence in the last equation 
will be uniform. 

Proof. The only part of the Corollary that does not follow directly from the preceding Theo- 
remEJis the uniformity of the convergence and that q will be harmonic on each of the intervals 
[ai,ai+i] for i < m. However, monotone convergence on compact sets preserves harmonicity 
and is always uniform (cf e g Meyer @ - or, more directly, Port and Stone 0 Theorem 3.9] 
if P is the Brownian semigroup). □ 

Lemma 7. In the preceding Corollary’s notation, q is the ininimal nonnegative fixed point 
of 1C. 

Proof. The proof partly copies the one for Lemma EJ Any nonnegative fixed point p of K, 
must be greater or equal g on [ao,a m ]. Therefore the monotonicity of 1C on [ao,a m ], implies 

Vn S No p = JC n p > IC n (g V 0) on [ao, a m ], 


yielding 


P > sup K n (g VO ) = q on [ao, a m ]. 

n£ No 


□ 
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3.4 Application 

Example 1 (Bermudan vanilla calls and puts on an asset that pays dividends at 
inflation rate in the Black-Scholes model). Assume 

P := {Pt)t >0 - = {ynt,a 2 t * ')t>0 ' 

where 



thus P can be perceived as the semigroup associated to the logarithmic price process under 
the risk-neutral measure in the one-dimensional Black-Scholes model). We will assume that 
(possibly after re-scaling the time scale) a = 1. Define 

go := —K + exp, g\ : K — exp 

(the payoffs on exercise of a one-dimensional call and a one-dimensional put option with 
strike price K, respectively). The infinitesimal generator of the Markov semigroup P is 

Let us, for the remainder of this Example, assume that r = 5. Then we have Lexp = 0 
and of course LI = 0 (1 denoting the constant function 1 : x i—» 1). Therefore, by Lemma [H 
the set of harmonic functions equals the space of all linear combinations of 1 and exp, which 
is the set of all real multiples of the exponential function with additive constants. Thus, both 
go and gi will be harmonic. 

Now suppose we are given support abscissas ao < ■ ■ ■ < a m as in the previous paragraphs. 

In order to obtain the setting of Theorem^ for g = go, set ho = a + exp for an arbitrary 
a > —K. Then ho will be a harmonic function dominating go = exp —K on R. Also, the 
constant function Co = minr ra0;O i (exp — K ) = e a ° —K will be harmonic and trivially less than 
or equal to go on [ao,a m ]. Then the conditions of Theorem |7| are satisfied for c = Co, g = go 
and h = ho, and the Theorem as well as its Corollary may be applied to establish the existence 
of a unique minimal fixed point for K. and its characterisation as the approximate price of 
the corresponding perpetual Bermudan option with exercise mesh t. - the approximation being 
based on piecewise harmonic interpolation as outlined in the introduction to this paper. In 
order for K. to be financially meaningful, however, we must choose a o > ln/\ to ensure that 
Co, the pointwise lower bound on the elements of Q, is nonnegative. 

Finally, in order to get the situation of Theorem |7] also for g = gi, we now define hi = 
K — exp to get a harmonic function dominating g\ = A'— exp on all of R. Also, the constant 
function ci = miiir ao , a i ( K — exp) = K — e“ m will be harmonic and trivially less than or 
equal to g i on [ao,a m ]. Then the conditions of Theorem^ are satisfied for the choices of 
c = Ci, g = g± and h = h\, and both the Theorem and its Corollary may be applied to 
see that again there is a minimal fixed point for /C coinciding with the approximate price 
of the corresponding perpetual Bermudan option with exercise mesh t. In order for K. to be 
financially meaningful, we must this time choose a m < In K to ensure that all elements of Q 
are pointwise nonnegative. 

4 Reduite-based approximation of Bermudan option 
prices 

Suppose P is a Markov semigroup on R d (d G N) and L is the infinitesimal generator of P. We 
will call a function / : R d —> R subharmonic if and only if 

Vt>0 P t f>f 

holds pointwise. A function / : R d —► R will be called superharmonic if and only if —/ is 
subharmonic, and / : R d —► R will be called harmomic if it is both super- and subharmonic. 


Let U denote the operator of upper-semicontinuous regularisation, that is, for all functions 

/ : R d -> R, 

Uf = inf {£ > f : £ : R d —■> R subharmonic} 

(of course, this is a priori only defined as a function taking values in RU{-oo}). Consider 

a harmonic function h : R d —> R and a closed (and therefore F a ) set B and define the reduite 
operator 1Z = lZh,B on the set of all subharmonic functions / : R d —> R dominated by h via 

7 Zf := U (sup {£ < h : £ : M d —> R subharmonic, l < / on B}) . 

It is a well-known result from potential theory (cf e g the work of Paul-Andre Meyer [5-. Theoreme 
T22]) that there will be a greatest subharmonic function dominated by / on B and that this 
function will be equal to IZf. Moreover, we have that / = 7 Zf on B except on a set of potential 
zero, in probabilistic/potential-theoretic jargon 

/ = IZf q.e. on B, 

where ‘‘q.e." is, as usual, short-hand for ‘‘quasi-everywhere”. Now define 

Q := {/ < h : f : R d —■> R subharmonic} . 

Then our definition of the reduite operator TZ implies IZf < h (as h is dominating the function 
whose upper-semicntinuous regularisation is, according to our definition, the reduite IZf of /) 
and our potential-theoretic characterisation of the reduite - as the greatest subharmonic function 
dominated by / on B - ensures the subharmonicity of IZf. Therefore, 

7 Z-.Q^Q. 

We also have that U is monotone (in the sense that for all /o < /i, Ufo < W/i) so that 1Z must 
be monotone as well (from the C-monotonicity of sup and the definition of 1Z). 

Hence 

Lemma 8. Adopting the notation of the preceding paragaph, 1Z : Q —> Q and whenever 

fo < fi, TZfo < TZfi- 

Let g : M. d —> R be a subharmonic function such that g < h and let r > 0. The next step is 
going to be the consideration of the following family of operators: 

<f> t ,f^e- rt P t f\lg 

for t > 0. If / < h , Ptf < Pth = h for all t > 0, since the operators Pt are positive and linear, 
and h was assumed to be harmonic. Thus, since g < h and r > 0, one must have <ptf < h 
for all f < h and t > 0. Moreover, the operators Pt preserve subharmonicity and the maximum 
of two subharmonic functions is subharmonic again, therefore 4 >tf must be subharmonic for all 
subharmonic /. Finally, since P t is monotone, </> t has to be monotone for all t > 0. Summarising 
this, we obtain 

Lemma 9. Using the notation introduced previously, 4t ■ Q Q an d whenever fo < f\, 
<t>tfo < <t>tfi for all t > 0. 

As a consequence, we derive from the two Lemmas0and O the following: 

Corollary 3. If we define /C* := TZo 0 t (adopting the notation of the previous paragraph), 
we have TC* : Q —> Q and whenever fo < fi, ICtfo < ICtfi. 

Corollary 4. The map f i—> Aft/V0 is a sound iterative Bermudan option pricing algorithm 
for the payoff function g V 0 (in the sense of Definition QJ). 

This already suffices to prove the following 
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Theorem 2. Let t > 0. Then for all n £ No 


/C t n+1 (gV0) >ICt n (gV0). 


( 1 ) 


Furthermore, 

q := sup tC t n {g V 0) 

n£ No 

(which a priori is only defined as a function with range in RU {+oo},) is an element of Q 
and indeed is the least nonnegative fixed point of ICt. 

Proof. 1. Relation 0 follows from the fact that IC t is a sound algorithm and Remark 0 

2. Since IC t maps Q to itself, the whole sequence (IC t n (g V 0)) n€No is bounded by h. This 
entails q < h as well. Applying Beppo Levi’s Theorem on swapping sup and / -dp 
- for bounded monotonely increasing sequences of measurable nonnegative functions 
and an arbitrary measure p - to the measures Pt(-,x), x £ and the sequence 
(ICt (ff v o)) ngNo , we can exploit the subharmonicity of the functions ICt n (g\/ 0), n € No, 
to deduce 

Vx G P t q(x) = sup P t ( K, t n (g V 0)) (x) 

n£ N 0 

> sup K.t n (g V 0)(x) = q(x), 
nG No 

which is the subharmonocity of q. As we have already seen, q < h, so q £ Q. 

3. If we employ Beppo Levi’s Theorem again, we can show that ICt and sup ngNo commute 
for bounded monotonely increasing sequences of functions. Thereby 

IC t q = sup lC t IC t n (g VO) = sup K. t n (g V 0) = q. 
n£No "SN 


4. That q is the least nonnegative fixed point is seen as in the proof of Lemma 0 Any 
nonnegative fixed point p of ICt must be greater or equal g V 0. Therefore by the 
monotonicity of sup and ICt, 

sup IC t n p > sup )C t n (g V 0) = q. 

nGNo 71 GNo 


□ 


Example 2 (Bermudan call option with equidistant exercise times in t ■ No on 
the weighted arithmetic average of a basket in a special Black-Scholes model). 

Let (Pi,..., fid) G [0,1] be a convex combination and for simplicity , assume that the assets 
in the basket are independent and each follow the Black-Scholes model with one and the 
same volatility a\ = •■■ = ad ='■ cr, and let r > 0 be the interest rate of the bond. We 
may assume that, possibly after a linear change of the time-scale, a = 1. Then (Pt.) t >o = 

is the semigroup of this Markov (even Levy) basket. Then one has 
t> 0 



L = 




■ V 


(cf e g Revuz and Yor's exposition m, and for 


d 

9 : x l-> & exp ” K 

2=1 


we obtain 

^ = E (it+ (’•-§)) «*«■><) 


10 


which is pointwise nonnegative if and only if 


> max ; ; e{ i 1 A 2 

T ~ 2 “ 2 d + 2d ' 

Hence, if r is sufficiently large, g is subharmonic and we can apply the theory developed 
earlier in this paper, in particular Theorem 03 


5 Soundness and convergence rate of perpetual Bermu¬ 
dan option pricing via cubature 

When Nicolas Victoir studied “asymmetric cubature formulae with few points" [IjQj for symmetric 
measures such as the Gaussian measure, the idea of (non-perpetual) Bermudan option pricing 
via cubature in the log-price space was born. In the following, we will discuss the soundness and 
convergence rate of this approach when used to price perpetual Bermudan options. 

Consider a convex combination (op£ [0, l] d (that is, Y^k=i a k = 1) and 
£ M. d . Then there is a canonical weighted arithmetic average operator A associ¬ 
ated with a, x given by 

m 

V/£R r Af = Y J *kf{-~x k ). 

fc =i 

Now suppose c £ (0,1), g , h : — > R, Ag > g, Ah = h and 0 V g < h. Define an operator V 

on the cone of nonnegative measurable functions by 

V : f ^ (c- Af)V g. 

We should stress that we explicitly allow for g to take negative values. 

Lemma 10. Adopting the previous paragraph’s notation and setting 

Q:={f<h : Af > / > 0} , 
we have the following properties of V and Q: 

1. T> is monotone (ie Vf < Vg whenever f < g). 

2. AT) — V is nonnegative on Q. 

3. V :Q ->Q. 

4- {g v o) £ Q. 

Proof. 1. Since A is positive and linear, thus monotone (in the sense that for all /o < /i, 
Afo < Afi), it follows that I? is a composition of monotone maps, thus monotone as 
well. 

2. Whenever Af > /, the monotonicity of A and our assumption Ag > g imply 

ADf > A(cAf) V Ag> ( cAf) V g = Vf. 

Hence AVf > Vf for all / £ Q. 

3. Because h is nonnegative and g V 0 < h and due to the monotonicity of V, we have for 
all nonnegative / < h, 

Vf < Vh = cAh V g = chV g = ch V0 V g <h, 

thus Vf < h for all / £ Q. Also, whenever Af > / > 0, one gets cAf > 0 and therefore 
in particular Vf > 0 for all / £ Q. However, we have already shown that AVf > Vf 
for all / £ Q. Summarising this, we arrive at Vf £ Q for every / £ Q. 
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4. Because of our assumption Ag > g and the monotonicity of A, we not only have 
A(g VO) > 0, but also A(gVO) > g, therefore A(gVO) > 5 VO > 0. However, by another 
assumption, g V 0 < h. Therefore (jVO)eQas claimed. 

□ 


Corollary 5. The map f 1 —> Vf VO is a sound iterative Bermudan option pricing algorithm 
for the payoff function g V 0 (in the sense of Definitional. 

Lemma ^21 suffices to prove 

Theorem 3. For all n £ No, 

V n+1 {gV0) >2?"(<?V0) =:q n . (2) 

Furthermore, 

q ■= lim D n {g V 0) = sup V n (g V 0) € Q 

n 1,00 n£N 0 

and q is the smallest nonnegative fixed point ofV. 

Proof. 1. Relation 0 follows from the soundness of D and Remark 0 

2. In Lemma HU we have not only seen that (g V 0) £ Q and V is closed under Q, but 
also that AVf > Df for all / £ Q. Hence 

VneNo A (D n (g V 0)) > V n (g V 0), 


and therefore 


Aq = sup A (V n {g V 0)) > sup V n {g V 0) = q, 

n£No nENo 

which means Aq > q > 0. 

Again because V maps Q itself and (j V 0) £ Q , we have that the whole sequence 
{D n (g V 0)) ngNo is bounded by h. This entails q < h as well. 

As we have already seen, Aq > q > 0, so q £ Q. 

3. Since A is a weighted arithmetic average operator, A and sup raeNo commute in the sense 
that A (sup n /„) = sup n Af n for increasing sequences of functions (/ ra ) n eN 0 - Hence, 
whenever (/ n )neN 0 is increasing, 

V ( sup fn ) = A ( sup fn ) V g = ( sup Af n ) V g = sup (Afn V g) = sup Df n , 

\ n / \ n / \ n J n n 

i e V and sup commute for bounded monotonely increasing sequences of functions. 
Thereby 

Vq = sup VV n {g V 0) = sup V n (g V 0) = q. 

raeNo n£N 

4. Just as in the proof of Lemma0 we see that q is the minimal nonnegative fixed point. 
For, any nonnegative fixed point p of V must be greater or equal g V 0. Thus, by the 
monotonicity of sup and D, 


sup V n p > sup V n (g V 0) = q. 

nENo nENo 


□ 


Lemma 11. Using the previous Theorem’s notation, we have for all x £ R d and n £ No, if 
q n +i(x) = g(x), then q n {x) = g{x). 

Proof. By the monotonicity of the sequence (o f n)neN 0 (Theorem 0, we have 

g{x) < qo(x) < q n (x) < q n+1 (x). 


□ 
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Theorem 4. For all n € N 


lkra+1 9n|lc°(R d ,R) — C ' ll^n Qn-1 ||<7 0 (R d ,R) ' 

Proof. The preceding Lemma, fill yields 

ll^n+l 9ra|lc , °(R d ,R) = II^JT+l ~ Qn llc 0 ({c-Aq„>g},R) 

= ||c • Aq n - ((C • Aq n _ i) V ff)llc°({c-A g „>g},R) 

via the definition of r/,; + i as ( cAqi ) V g for i = n and i = n + 1. But the last equality implies 

lkn+1 — 9" llc°(R d ,R) — ll c ' Aq n — C ■ Aq n _i |lc'0({ c - J 4q„ >g },R) 

— II c ■ Aq n — c ■ Aq n _i ■ 

Since A is linear as well as an T°°-contraction (and therefore a C°-contraction, too), we 
finally obtain 

||<Zra+l — Qn llc°(R d ,R) — c II A (q n ~ 9n —1) || C'°(R d ,R) — c 119™ — Qn — 1 llc°(R d ,R) ' 

□ 

Example 3 (Bermudan put option with equidistant exercise times in t ■ No on 
the weighted arithmetic average of a basket in a discrete Markov model with a 
discount factor c = e~ rt for r > 0 ). Let / 3 i,..., f3d € [ 0 , 1 ] be a convex combination and 
assume that A is such that 


V*€{l,...,d} J2 ake ~ (Xk)i = 1> | 3 ) 

k =1 

then the functions 

d 

g : x i->- K - ^2 Pi ex P ( x i) 

2=1 

and h := K (where K >0) satisfy the equations Ah = h and Ag = g, respectively. Moreover, 
by definition g < h, thus 0 V g < h. Then we know that the (perpetual) Bermudan option 
pricing algorithm that iteratively applies V to the payoff function g VO on the log -price space, 
will increase monotonely and will have a limit which is the smallest nonnegative fixed point 
ofV. Moreover, the convergence is linear and the contraction rate can be bounded by c. 

The condition m can be achieved by a change of the time scale (which ultimately leads 
to different cubature points for the distribution of the asset price) 

One might also be interested in determining the convergence rate for the approximation of non- 
perpetual American option prices based on non-perpetual Bermudan option pricing via cubature. 
After proving a series of Lemmas we will end up with a Theorem that asserts linear convergence 
and also provides bounds for the convergence factor. 

From now on, c and A will no longer be fixed but their role will be played by e~ rt and P* (for 
t e sNo where s > 0 shall be fixed) respectively, where r > 0 and (P s . m ) mgNo describes a Markov 
chain on (By the Chapman-Komogorov equation this is tantamount to Vs,i > 0 P s Pt = 
Ps+t). 
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